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We consider a heavy meson moving slowly through high temperature non-abelian 
plasmas. Using a simple dipole effective Lagrangian, we calculate the in-medium 
mass shift and the drag coefficient of the meson in N = 4 Super Yang Mills theory 
at weak and strong coupling. As anticipated, in the large N limit the mass shift is 
finite while the drag is suppressed by 1/Y 2 . After comparing results to perturbative 
QCD estimates (which are also calculated), we reach the conclusion that relative to 
weak coupling expectations the effect of strong coupling is to reduce the momentum 
diffusion rate and to increase the relaxation time by up to a factor of four. 



I. INTRODUCTION 

The energy loss of heavy quarks and quarkonia in media has been a subject of intense 
experimental interest [1-7]. The suppression of charm and bottom quarks observed at RHIC 
motivated several groups to utilize the guage-gravity duality [8-11] to compute the drag of 
fundamental heavy quarks in Af = 4 Super Yang-Mills at strong coupling [12-14]. The goal 
of this paper is to extend these strong coupling calculations of drag and diffusion from heavy 
quarks to heavy mesons. 

The motivation for this effort is twofold. First, future experiments at RHIC promise to 
measure the elliptic flow of J/ip mesons, and it is important to support this experimental 
program with theoretical work. To this end, various groups have studied the thermal proper- 
ties of heavy mesons within the context of the AdS/CFT correspondence [15-18]. However, 
in spite of this progress, the transport properties of these mesonic excitations are not well 
understood. Although the kinetics derived in this paper are not directly applicable to the 
heavy ion experiments, we believe that the results do hold some important information for 
phenomenology. 



*Electronic address: kdusling@gmail.com 
'''Electronic address: jke@mppmu.mpg.de 
Electronic address: kaminski@mppmu.mpg.de 
§ Electronic address: rust@mppmu.mpg.de 
^Electronic address: derek.teaney@stonybrook.edu 
"Electronic address: young@tonic.physics.sunysb.edu 



2 



The second motivation for this work is theoretical. After the quark drag was computed 
using the correspondence, it was realized that the drag of quarkonia is zero in a large N c 
limit [15, 16, 18]. In the high temperature phase (when the temperature is comparable 
to the mass of lowest meson state) the vector spectral function has been computed and it 
shows a rich dynamical picture of meson melting [19-23]. However, in the low temperature 
phase describing heavy mesons, the spectral function is usually described by a sequence of 
states with zero width. Since within a thermal environment the drag and diffusion of these 
mesonic states is certainly not zero, it remained as a theoretical challenge to compute the 
kinetics of these states using the AdS/CFT setup. 

Rising to this challenge, the width of AdS/CFT mesons due to scattering with surround- 
ing heavy quarks (or anti-quarks) was recently determined by extending the analysis of 
meson melting to finite baryon density [19] and by studying string worldsheet instantons at 
zero baryon density [24]. In general the meson width determined in this way is suppressed 
by the density of heavy quarks. (At zero baryon density the width is suppressed by the 
thermal population of heavy quarks.) In contrast, we are concerned with the thermal width 
which is finite at zero density and infinite A and captures the rescattering between the meson 
and the surrounding M = 4 medium. 

This work will focus on heavy mesons where the binding energy is much greater than 
the temperature. While in perturbation theory the constraint on the binding energy reads 
rriqV 2 3> T, in the strongly coupled Af = 4 theory the constraint is 3> T. In this tight 
binding regime, mesons survive well above T c and the meson width is sufficiently narrow to 
speak sensibly about drag and momentum diffusion. 

For real charmonium (bottomonium) the binding energy can be estimated from mass 
splitting between the 2s and Is (3s and Is) states, AM 2 ^ ls ~ 589 MeV and AM^ s _ ls ~ 
895 MeV respectively [25] . Therefore it is not really clear that real quarkonia above T c ~ 
170 <-> 190 MeV [26, 27] can be modeled as a simple dipole which lives long enough to 
be considered a quasi-particle. Indeed weak coupling hot QCD calculations of the spectral 
function show that over the temperature range g 2 M <-> gM, the concept of a meson quasi- 
particle slowly transforms from being well defined to being increasingly vague [28-32] . There 
is lattice evidence based on the maximal entropy method (which is not without uncertainty) 
that J/ip and T survive to 1.6 T c and ~3T C respectively [33-38]. However, model potential 
calculations which fit all the Euclidean lattice correlators indicate that the the J/ip and T 
survive only up to at most 1.2 T c and 2.0 T c respectively [39, 40]. Clearly, the word "survive" 
in this context is qualitative and means that there is a discernible peak in the spectral 
function. Given these facts, the assessment of the authors is that the dipole approximation 
might be reasonable for T(IS') but poor for charmonium states and other bottomonium 
states. 

Within the context of guage gravity duality, mesons are studied by exploiting convenient 
generalizations of the AdS/CFT correspondence in which fundamental flavor degrees of 
freedom are added by inserting additional probe branes into the geometry. In this paper we 
will use the approach due to Karch and Katz [41] in which the additional quark flavors are 
obtained by adding D7 brane probes to the original D3 brane setup. On the gravity side, 
the probe branes wrap a subspace which asymptotically near the boundary is AdS 5 x S* 3 . 
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The meson spectrum for the dual M = 2 supersymmetric gauge theory was first calculated 
in [42] by considering fluctuations of the D7 branes embedded. Restricting this result to 
fluctuations with vanishing angular momentum on the S 3 , the meson spectrum is given by 

M=^2v/(n + l)(n + 2), (1.1) 

with m q the quark mass determined by the separation between the D3 and D7 branes, A the 
't Hooft coupling and n the radial excitation number. Subsequently, a gravity dual of chiral 
symmetry breaking has been obtained in [43] by embedding a D7 brane probe into a deformed 
non-supersymmetric gravity background with non-trivial dilaton [44] . In this case there is a 
Goldstone boson in the meson spectrum. The thermodynamics of mesons has been studied 
within gauge/gravity duality by embedding a D7 brane probe into the AdS-Schwarzschild 
black hole background. Within the deconfined phase, there is a new fundamental first order 
phase transition which corresponds to meson melting [21, 43, 45-47]. For a review on mesons 
in the AdS / CFT correspondence see [48] . 

In this paper we consider a heavy meson moving slowly through the medium. We per- 
form both a perturbative QCD and a strong coupling M = 4 SYM computation. For both 
approaches we first calculate the in medium meson mass shift, which determines the polariz- 
abilities of the meson. As expected from the dipole effective theory, the mass shift scales as 
T 4 /Ag, with T the temperature and A B the inverse size of the meson. In the perturbative 
calculation, A B is the inverse Bohr radius, while in the AdS/CFT computation the meson 
mass plays this role. In the M = 4 field theory the dipole effective Lagrangian couples the 
heavy meson to the stress tensor and the square of the field strength F 2. In AdS/CFT we 
obtain these couplings from the linear response of the meson mass to switching on a black 
hole background or a non-trivial dilaton flow, respectively. For the dilaton flow we consider 
the D3 + D(— 1) gravity background of Liu and Tseytlin [49]. This background and the 
AdS-Schwarzschild background allow for an analytic calculation of the meson polarizabili- 
ties. 

Using these polarizabilities we subsequently calculate the momentum broadening k and 
the drag coefficient i] D . This requires the calculation of two-point functions involving gradi- 
ents of the stress tensor and the field strength squared. Within gauge/gravity duality, these 
are obtained by considering graviton and dilaton propagation through the AdS-Schwarzschild 
black hole background. 

An outline of the paper is as follows. First, in Section II we review the computation 
of drag and diffusion of heavy QQ bound states within the setup of pertubative QCD. 
This will outline a two step procedure to determine the drag coefficient at strong coupling. 
The first step is to determine the in medium mass shift (which is finite at large N c ) which 
determines the polarizabilities of the meson. This is done in section Section III. The second 
step is to compute the force-force correlator on the meson using the previously computed 
polarizabilities. This determines the drag and diffusion coefficient as reviewed in Section IV. 
Finally we compare our results to perturbation theory and reach some conclusions for the 
RHIC experiments in Section V. 
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II. DIFFUSION OF HEAVY MESONS IN PERTURB ATIVE LARGE N FIELD 

THEORIES 

A. Diffusion of Heavy Mesons in Perturbative large N QCD 

The interactions of a heavy meson with the QCD medium is well described by a dipole 
approximation. This physical approximation has been formalized in the language of heavy 
meson effective Lagrangrians which we will adopt [50]. This perturbative scheme relies on 
the large mass of the meson relative to the external momenta of the gauge fields but does 
not rely on the smallness of the coupling constant. It was used previously to make a good 
estimate for the binding of J/ip to nuclei [50]. 

The heavy meson field <p v describes a (scalar) meson which has a fixed velocity v M = 
(7,7V). Then the effective Lagrangian for this meson field interacting with the gauge fields 
is 

C E , , C B 

^cflf = ~(pvlV ■ 0(p v + 

with 



£ cff = -4>\iv ■ d<p v + ^Oe^ + ^IOb^ , (2.1) 



O e = - l -G^ A G a vA v^ v , O b = ^G a ? A G af3 A - \p mA Ga A v ll v v . (2.2) 

G^ v is the non-Abelian field strength of QCD, and ce and cb are matching coefficients (po- 
larizabilities) to be determined from the QCD dynamics of the heavy QQ pair. In inserting 
a factor of 1/N 2 into the effective Lagrangrian we have anticipated that the couplings of the 
heavy meson to the field strengths are suppressed by N 2 in the large N limit. 

In the rest frame of a heavy quark bound state, v = (1, 0, 0, 0), the operators Oe and Ob 

are 

O e =^E a -E a , Ob = \b a -B a , (2.3) 

where E A and B A are the color electric and magnetic fields. If the constituents of the dipole 
are non-relativistic it is expected that the magnetic polarizability cb is 0(v 2 ) relative to the 
electric polarizability. For heavy quarks (where cb is neglected) and large N these matching 
coefficients were computed by Peskin [51, 52] 

ce = H| , c B = . (2.4) 

Here A B = l/a Q = (m g /2)Ci?o; s is the inverse Bohr radius for a QQ bound state. It is finite 
at large N since with Cf — N/2 and finite A we have A B = m q \/l6ir. We will assume that 
cb = is zero in this section and subsequently generalize our results to Af = 4 theory. 

This effective Lagrangian can be used to compute both the thermodynamics and kinetics 
of the heavy meson state. First, to evaluate the in medium mass shift one simply uses first 
order perturbation theory 5M = (Hi) = — (£1), yielding 

A J 45 • 
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In the second line we have simply calculated the expectation value (Oe) t = ^N 2 T 4 in a 
free gluon gas and used Eq. (2.4). 

The importance of this result is that it is finite at large N and that it is in general 
suppressed by (T/A^) 3 , i.e. suppressed by powers of the hadron scale to the temperature. If 
higher dimension operators were added to the effective Lagrangian their contributions would 
be suppressed by additional powers of T/A B . At strong coupling, we will use the AdS/CFT 
correspondence to determine the M — 4 polarizabilities from the mass shift. 

We turn next to the kinetics of a heavy dipole in the medium. For timescales which are 
long compared to medium correlations, we expect that the kinetics of the heavy meson is 
described by Langevin equations 

^ = m - voPi , mMf)) = ^5{t - a) . (2.6) 

Here £j is a random force with second moment k and i] D is the drag coefficient. The drag 
and fluctuation coefficients are related through the Einstein relation 

^ = 2Wr- (2 ' 7) 

The Langevin equation is valid for times which are long compared to the inverse temperature 
but short compared to the lifetime of the quasi-particle state. 

We can use the effective dipole Lagrangian to calibrate the noise of the stochastic evolu- 
tion, k. The microscopic equations of motion for a heavy particle in the medium are 

^T = ^W, (2-8) 

where T is a phenomenological force to be specified below. We then compare the response 
of the Langevin process to the microscopic theory. Over a time which is long compared to 
medium correlations but short compared to the time scale of equilibration we can neglect 
the drag and equate the stochastic process to the microscopic theory 

Jdt J dt' (&(*)&(*')) = (time) x /t<% = J dt J dt' (^(0^(0) ■ ( 2 - 9 ) 

In a rotationally invariant medium we have 

K=\fto (F(t)F(0)) . (2.10) 

In the present context we identify the force with the negative of the gradient of the interaction 
Hamiltonian Hi = —Ci , 

F{t) = J d 3 x0t(x,t) [c B VO £ (x,t)]^(x,i) , (2.11) 

which is the usual form of a dipole force averaged over the wave function of the meson. 
The fluctuation dissipation theorem relates the correlation function in Eq. (2.10) (with the 
specified time order of operators) to the imaginary part of the retarded force-force correlator 

1 — 2T 

k = - lim lmG R (u) , (2.12) 

3 ^-+0 uj 
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FIG. 1: Dipole scattering graph which causes drag and diffusion of heavy mesons in QCD and 
N = 4 SYM. The QCD correlator which encodes this physics is given by Eq. (2.14). 

where the retarded correlator is 

Gr — —i Jdte +iut 9(t)([r(t),r(0)]) . (2.13) 

Integrating out the heavy meson field as discussed in detail in Ref. [13], which treated the 
heavy quark case, we obtain a formula for the momentum diffusion coefficient 



-^imcri-.q) 



(2.14) 



3N* J (27r)3 q 
with the retarded OeOe correlator given by 

G% E ° E (iu,q) = -iJ d 4 xe + ^-^9(t) ([0 E (x,t) , O E (0,0)]) ■ ( 2 -15) 

We can understand this result with simple kinetic theory. Examining the Langevin dy- 
namics we see that 3k is the mean squared momentum transfer to the meson per unit time. 
The factor of three arises from the number of spatial dimensions. In perturbation theory 
this rate is easily computed by weighting the transition rate for any gluon in the bath to 
scatter with the heavy quark by the square of the momentum transfer, 

3K = / A to, ^ + v) i 2 - p + p-> ■ 

Here p is the incoming gluon, p' is the outgoing gluon and q is the momentum transfer 
q = p — p' as indicated in Fig. 1. \Ai\ 2 is the gluon meson scattering amplitude computed 
with the effective Lagrangian in Eq. (2.1) and summed over colors and helicities of the 
incoming and outgoing gluon 

\M\ 2 = ^-y (l + cos 2 (V)) • (2-17) 

Alternatively (as detailed in Appendix A) we can simply evaluate the imaginary part of the 
retarded amplitude written in Eq. (2.14) to obtain the same result. 
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For QCD the integrals written in Eq. (2.16) are straightforward and yield the following 
result for the rate of momentum broadening 

1 2 64tt 5 9 

T 3 /7rT\ 6 501767r (2 ' 18) 



N 2 \A B J 1215 ' 

The high power of temperature T 9 arises since the dipole cross section rises as uj a . The 
matching coefficient ce is directly related to the mass shift of the dipole and encodes the 
coupling of the long distance gluonic fields to the dipole. By taking the ratio between the 
momentum broadening and the mass shift squared we find a physical quantity which is 
independent of this coupling 

" lT128 ° (2.19) 



(5M) 2 N 2 3 

The large numerical factor 1280/3 originates from the cross section which grows as uj a . A 
similarly large factor appears in Af = 4 SYM as discussed below in greater detail. 



B. Linear perturbations of N = 4 Super Yang- Mills theory 

Our aim is to calculate the heavy meson diffusion coefficient k from gauge/gravity du- 
ality. This requires the calculation of the two-point correlators as well as of the associated 
polarizabilities in Af = 4 Super Yang-Mills theory. 

The same formalism used in the preceding section can be used for Af = 4 SU (N) Super 
Yang-Mills theory. In general all operators in Af = 4 which are scalars under under Lorentz 
transformations and SU (4) R-charge rotations will couple to the meson at some order. The 
contribution of higher dimensional operators is suppressed by powers of the temperature to 
the inverse size of the meson. The lowest dimension operator which could couple to the 
heavy meson field is O x ^ = trX' l X l , where X 1 denotes the scalar fields of Af = 4 theory. 
However the anomalous dimension of this operator is not protected, and the prediction of 
the supergravity description of Af = 4 SYM is that these operators decouple in a strong 
coupling limit [11]. The lowest dimension gauge invariant local operators which are singlets 
under SU (4) and which have protected anomalous dimension are the stress tensor T^ u which 
couples to the graviton, minus the Lagrangian Opi = —Cn=4, which couples to the dilaton 1 
and the operator O ff = tr F^F^-^. . ., which couples to the axion. An interaction involving 
F p breaks CP which is a symmetry of the Lagrangian Af = 2 hypermultiplet of the Af = 4 
SYM gauge theory. Thus interactions involving O ff can be neglected. 

Summarizing the preceding discussion, we find that the effective Langrangian describing 
the interactions of a heavy meson coupling to the operators in the field theory is 

£ cff = -<f> v (x,t)ivd<l> v (x,t) + ^0i(x,t)T^^0„(x,t) + ^<t>\{x,t)0 F 2<t> v {x,t) , (2.20) 



Since we can add a total derivative to the Lagrangian, the operator — L is ambiguous. The precise form 
of the operator coupling to the dilaton is given in Ref. [53]. We neglect this ambiguity here. 
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which is a linear perturbation of Af = 4 Super Yang-Mills theory by two composite operators. 

The polarization coefficients ct, c f will be determined below from meson mass shifts in 
gauge/gravity duality. This requires breaking some of the supersymmetry. For the contri- 
bution of the energy-momentum tensor, this is achieved by switching on the temperature. 
Then, the mass shift of the meson is given by expectation value of the stress tensor 

5M = -^<T 00 > , (2.21) 

In gauge/gravity duality this is achieved by considering the AdS-Schwarzschild black hole 
background where (Opt) = 0. In contrast, for the meson response to (Of 2 ) we consider a 
background self-dual gauge configuration where (0 F 2) ^ while (T^ u ) = 0. As can be seen 
from the supersymmetry transformations of the Af = 4 fermions, such a background breaks 
the supersymmetry to Af = 2. The mass shift of a heavy meson is then 

5M = -^ 2 {O f2 ) , (2.22) 

where (Of 2 ) = ( _ £) is the expectation value of the Af = 4 Lagrangian in the Af = 2 sym- 
metric background configuration. A suitable dilaton solution corresponding to this selfdual 
configuration been given by Liu and Tseytlin [49] and will be used below to determine c F . 

As explained above for the QCD case, to determine the kinetics of the heavy meson in 
the Af = 4 background we identify the force on the heavy meson as minus the gradient of the 
interaction Hamiltonian Hj = —Cj averaged over the meson wave function. With Eq. (2.20) 
we have 

F(t) = J d 3 x0t(x,t) V \^T^u,u v + ^O f2 ] 0,(x,t) . (2.23) 

Then integrating out the heavy meson fields as above, the force- force correlator in Eq. (2.12) 
becomes 



f d 3 q q 2r/ CT x2_2T TT ,^(CfV-2T p2p2 



(2.24) 



where the retarded correlators are 



G T R T = -i Jd 4 xe +iult - l ^6(t)([T 00 (x,t),T 00 (0,0)]) , (2.25) 
G F R F = -i Jd 4 xe +iu)t ~^6(t)([O F 2(^t),O F 2(0,0)}} . (2.26) 

In writing Eq. (2.24) we have implicitly assumed that there is no cross term between F 2 
and T^v^Vy. In the gauge/gravity duality this is reflected in the fact that at tree level in 

supergravity, ~ ffg)^) = °' 

In summary, we first will determine the polarizabilites c F and ct from the mass shifts of 
the meson in two different backgrounds using Eq. (2.21) and Eq. (2.22). Subsequently we 
will compute the correlators in Eq. (2.25) for Af = 4 theory at finite temperature. Finally we 
will put the results together using Eq. (2.24) to deduce the rate of momentum broadening. 
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III. DETERMINING MATCHING COEFFICIENTS WITH MASS SHIFTS 
A. Backgrounds dual to finite temperature and field strength 

1. Finite temperature background 

The gravity background dual to M = 4 SYM theory at finite temperature is given by 
the AdS-Schwarzschild black hole with Lorentzian signature (see e.g. [54]). This back- 
ground is needed below both for calculating the necessary two-point correlators (T 00 T 00 ) 
and (Of^Of^), as well as for obtaining the meson polarizability, ct- 

We make use of the coordinates of [43] to write the AdS-Schwarzschild background with 
Lorentzian signature as 

ds 2 = ^ (- 1 dt 2 + /dx 2 ) + ^ (d, 2 + g'dQt + dwl + dwl) , (3-1) 

with the metric df2 2 of the unit 3-sphere, and 

r 'H = TttR 2 , R A = Ang s Nt s , A — ing s N, g 2 YM = g s . 
This spacetime has a horizon at % which is determined by th as 

w H = (3.3) 

and the boundary is reached at asymptotically large w. With r H = we obtain AdS 5 x S 5 . 
In section IV we will work in a coordinate system with inverted radial AdS coordinate u 
used in e.g. [54]. In these coordinates, the metric reads 

ds 2 = (-/(«) dt 2 + dx 2 ) + du 2 + R 2 dn 2 5 , (3.4) 

with f(u) = 1 — u 2 and df2 2 the unit 5-sphere metric. 



2. Dilaton background 



A non-trivial dilaton background that is dual to a field configuration with {Of 2 ) 7^ 
and {T^ v ) = has been given by Liu and Tseytlin [49] and consists of a configuration 
of D3 branes with homogeneously distributed D(-l) instantons. The type IIB action in the 
Einstein frame for the dilaton <£>, the axion C, and the self-dual gauge field strength F 5 = *F 5 
reads 



Sr 



—2~ [ dWx 

2k 10 J 



K 



„2$ 



{dCf 



1 

4^5! 



(F 5 ) 2 + 



(3.5) 
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The ten-dimensional Newtons constant is 

1 1 N 2 



(3.6) 



Solving the equations of motion derived from (3.5), Liu and Tseytlin [49] obtain the metric 

2 



^)VdaW+(^) (dr 2 + r 2 dfi 2 ) 



(3.7) 



d „2 _ «/2 J 2 _ */2 

u,5 string ~~ c us Emstem ~~ c 

and axion-dilaton solution 2 

- 1 + \ > C = -i (e"* - 1) . (3.8) 



e 



The expectation value (Of 2 ) is given by 



5S m N 



2 



B. Computing the polarization coefficients from meson mass shifts 

Heavy mesons are identified with fluctuations <p of a D7 brane embedded into the back- 
ground dual to the field theory under consideration. Stable embeddings are obtained if the 
D7 brane spans all Minkowski directions as well as the radial AdS coordinate and a 3-sphere 
in the remaining polar directions. Consider the metric (3.1) as an example. The D7 brane 
spans all directions except w$ and wq. The meson mass M is then obtained by solving the 
equation of motion for the fluctuations <p [42] . Read as an eigenvalue equation, the equation 
of motion gives the meson mass as the eigenvalues M to the corresponding eigenfunctions 
(p. The discrete values of M describe the Kaluza-Klein mass spectrum of mesons for any 
given quark mass. 

To see how this works, we outline this procedure for the vacuum case (T 00 ) = (Op*) = 
0, for which the meson spectrum was originally calculated in [42]. Subsequently we will 
introduce a non-zero (0 F 2) and (T 00 ), respectively. 

In the case of a D7 brane embedded in a ten dimensional background, the brane embed- 
ding is described by the location in the two directions transverse to the brane. We call these 
directions u>5 and wq. In general these locations depend on all eight coordinates of the 
eight dimensional D7 brane worldvolume and are determined by extremizing the DBI-action 

d 8 £e-*V^feth, h ah = —^G%, (3.10) 

where T7 is the D7 brane tension and G st is the string frame metric of the ten dimensional 
background with coordinates X^. It is related to the Einstein metric as in Eq. (3.7). The dis- 
tinction between the Einstein and string frame is ultimately important below. The pullback 



2 For the conventions note that in our notation q = ^qur, where q^x is used in the paper of Liu and 
Tseytlin [49]. 
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h contains the functions tus(^) and wq(£), which are determined by solving their equations 
of motion, derived from Sdbi- 

The background AdS^ x S 5 dual to (T 00 ) = (Of 2 ) — is obtained e.g. from (3.1) with 
th = 0. It is well known that for this background a probe brane embedding is given by the 
functions 

w 5 = 0, (3.11) 
w 6 = L, (3.12) 

where w§ and wq are the coordinates given in (3.1) and the constant L determines the quark 
mass m q = L/{2ti£ 2 ). Now we allow for small fluctuations <p 5 and y5 6 around this solution, 

w 5 = 2ir£ 2 s <p 5 (x, g) , (3.13) 
w 6 = L + 2tt£ 2 s <p 6 (x , g) . (3.14) 

By the symmetries of the setup, the fluctuations only depend on the Minkowski directions x 
and on the coordinate g, denoting the radial coordinate on the part of the D7 brane which 
is transverse to the Minkowski directions. The resulting equations of motion are identical 
for <f 5 and <p 6 . Using <p to denote any one of them, it was shown in [42] that a solution to 
the equations of motion may be obtained by separating variables with the ansatz 

V = V (g)e rM y\S*) , (3.15) 

where y e (S 3 ) are the scalar spherical harmonics on the S 3 wrapped by the probe D7 brane 

— * 

and k denotes a four vector. The resulting equation of motion for the radial function (p(g) 
may be obtained from (3.10). For £ = it reads as 

-d pP % V (p) = M 2 {p f +l)2 ip{p) ■ (3.16) 

Here we introduced the following dimensionless quantities 

P=f, M= T M, W = ^f, (3.17) 

and identified the meson mass squared M 2 with the square of the momentum four-vector k 
of the fluctuations, 

M 2 = -k 2 . (3.18) 

The eigenfunctions y? n solving the Sturm-Liouville equation (3.16) are given in terms of 
the standard hypergeometric function 2 F\-, 

MP) = (p2 + C " 1)n+1 2F 1 (-(n + 1); -n; 2; -p 2 ) , (3.19) 
where c n is a normalization constant such that 

no 

p 3 

d P ( p2 + iy Vn{p) <Pm(p) = $nm ■ (3.20) 
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The lowest mode ipo is given by 

'12 



<Po(p) = dnrr • (3-21) 



p 2 + l 

The corresponding eingenvalues M n to the functions ip n are given by 



M n = 2y/(n + l)(n + 2) . (3.22) 
We note that the mass of the lowest state with n = is 

M = -^2v / 2=^2v / 2, (3.23) 

which will appear frequently below. For a more detailed derivation of these results the reader 
is referred to [42]. 



1. Mass shift in the dilaton background 

Let us now calculate the polarizability c F which determines the change SM of the meson 
mass at a given value of the gauge condensate (Opt) with respect to the meson mass at 

5M = -^- 2 (O f2 ) . (3.24) 

To find cf we will determine the mass shift 5M and identify c F with the proportionality 
constant in front of (C^f 2 )- 

We are interested in the eigenvalues of fluctuations in the case of q oc (Cf 2 ) 7^ 0. The 
ten-dimensional background geometry dual to this scenario is given in (3.7) and the equation 
of motion for D7 brane fluctuations analog to (3.16) was derived in [55] to be 

-BSBrfp) = - v + 1)( / + V + 1)2) 9Mp) ' (3 ' 25) 

with the dimensionless q 

Q = J* ■ (3-26) 

To obtain analytical results, we consider the case of small q and linearize in this parameter. 
Therefore the equation of motion to solve is 

-d p p 3 dMp) = M 2 {p /^ l)2 <p(p) + A(pMp) , (3.27) 

where the operator A(p) is given by 

A (p) = - 4 vttf dp - (3,28) 



It is this term that describes the difference between the equation of motion at non vanishing 
background perturbation to (3.16), which is valid for q — 0. 
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To find the solution <po(p) corresponding to the lightest meson with n = we set up a 
perturbative expansion. Any deviation Sipo from the solution tpo of the case q = may be 
written as a linear combination of the functions ip n , which are a basis of the function space 
of all solutions, 

oo 

4>{p) = MP) + ^MP) , On < 1 , (3-29) 

n=0 

M 2 = Ml + 5Ml , 5Ml < 1 . (3.30) 



Plug this ansatz into the equation of motion (3.25), make use of (3.16) and keep terms up 
to linear order in the small parameters a n , q and 5Mq to get 



Q OO q q OO 



( p2 + 1)2 J2 a nMlVn(p) = 6M$ + V^o(p) + M Q 2 ^ °^n(p) + A(p)v9 (p) ■ 

(3.31) 

We now multiply this equation by <po(p)j integrate over p e [0, oo] and make use of (3.20) 
and (3.21) to see that 

OO 

5Ml = - [ dp MP) &(p)Mp) 

o (3-32) 



From 8M 2 = 2M SM we obtain 



where in the last step we used (3.23) for the mass and (3.9) and (3.26) to relate q and Op2. 
By comparison with (3.24) we identify the polarizability 

8 /2tt x ; 



c ' = sUJ ' (3 ' 34) 

2. Mass shift in the finite temperature background 

The calculation of the polarizability ct is completely analogous. We are now looking 
for the proportionality constant of meson mass shifts with respect to deviations from zero 
temperature, 

5M = --^<T°°> . (3.35) 

The background dual to the finite temperature field theory is the AdS black hole background 
given in (3.2). Notice that the black hole radius r H is related to the expectation value (T 00 ) 
by [56] 

(T 00 ) = -n 2 N 2 T 4 , r H = tiTR 2 . (3.36) 
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; r H = l \ 
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0.0 0.5 1.0 1.5 2.0 

P 

FIG. 2: Some embedding functions we(p) for different quark masses m q = linip^oo w§/ '(2tt£ 2 ). For 
a detailed discussion see e.g. [43]. 



Again we calculate the meson mass spectrum to identify the polarizability by comparison 
with (3.35). The embedding functions w^g) and wq(q) in this background are given by 

w 5 = , (3.37) 
w& = w 6 (g) , (3.38) 

where the quark mass is determined by m q = lim^oo wq/{2-k1 2 s ). The function Wq(q) has to 
be computed numerically [43]. Some examples of such embeddings are shown in Fig. 2. We 
introduce small fluctuations (p(p)e lkx in the w^, direction, 



w 5 -> w 5 (g,x) 



ikx 



(3.39) 



The linearized equation of motion for the fluctuations <p(g) in the limit of vanishing spatial 
momentum and M 2 = —k 2 can be derived from the DBI action (3.10) to be 







d e <p(g) 



1 + (d e w 6 y 



Q 



2{ Q 2 + wlf 



r%<P(Q) 



where we abbreviated 



(3.40) 



H 



{{q 2 + wD 2 



G = e 3 i 



Hi 



H 



16 (p 2 + w 2 ) 4 



(3.41) 



In the regime of small temperatures, we may linearize in rjj which is the leading order 
in th- Furthermore, as may be seen from figure 2, in the regime of a small temperature 
T compared to the quark mass m q , or respectively small ratios of r#/ lmip^oo wq(p), the 
embeddings become more and more constant. So for constant embeddings wq = L and up 
to order T 4 oc rjj the equation of motion simplifies to 



-d p p 3 d pV (p) = M 1 



(P + l) 



Mp) + A(p)<p(p) 



(3.42) 
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where we made use of the dimensionless quantities (3.26) and identify 



4L 4 (p 2 + l) 4 

For the lightest meson, the ansatz (3.29) this time leads to 



A(P) = V-kj^V^ 2 ■ (3-43) 



oo 



5Ml = - Jdpip {p)A{p)ip (p) 



o 

9 rjjMl 
40 L 4 ' 

Reinstating units and solving for SM leads to 

12 /2tt\ 3 1 



(3.44) 



From this we can read off the polarizability c F as 



12 / 2tt x 3 



*=5?UJ ■ < 3 ' 46 > 

IV. FINITE TEMPERATURE CORRELATORS 

According to Eq. (2.24) we need to compute the following correlators at finite temperature 

G T R T (u,q) = -i j d 4 x e +i "'- iq - x 9 (t)<[T 00 (x,t),T 00 (0,0)]) , (4.1) 

G F R F (u,q) = -i Jd*xe + ^-^9(t)([O F 2(x,t),O F 2(0,0)]) . (4.2) 

The calculational procedure for these two correlators is standard and has been discussed in 
[57, 58]. T 00 correlators are associated with graviton propagation and F 2 correlators are 
associated with dilaton propagation. 

On the gravity side both field correlators are computed in the black hole background (3.1) 
placing the dual gauge theory operator correlation functions at finite temperature. For 
simplicity in this section we work in the conventions and coordinates of [54]. We apply 
the method developed in [57, 58] and first applied in [54], in order to find the two-point 
Minkowski correlators as 

G R (u, q) = A(u) f(u, -k) d u f(u, k) . (4.3) 

The function f(u, k) relates the boundary and bulk values of a gravity field to each other. 
For example the dilaton field <3> is related to its value at the boundary bdy by 



<S>(u,k) = f(u,k)<p hd ^k), (4.4) 
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and is normalized to one at the boundary f(0,k) = 1. For metric fluctuations, &(u, k) is 
replaced by h 00 (u, k). The factor A(u) can be read off from the classical supergravity action 

1 f dud 4 xA(u) (<9 U $) 2 



Sr. 



+ 



The classical gravity action for the graviton and dilaton is obtained from (3.5) as 
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du d 4 x\/—g5 



(TZ-2A) - hd$) 2 + ... 



where 



N 2 



K 



10 



An 2 R 3 



So comparing to (4.5) we get 



2«§ 



The equation of motion derived from (4.6) in momentum space reads 



-<£>' + 



$ = , 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



uf(u) uf(u) 2 

with the function f(u) = 1 — u 2 , the dimensionless frequency to = uj/2ttT and spatial mo- 
mentum component q = qj2-nT . The equation of motion (4.9) has to be solved numerically 
with incoming wave boundary condition at the black hole horizon. Computing the indices 
and expansion coefficients near the boundary as done in [59, 60], we obtain the asymptotic 
behavior as linear combination of two solutions 

$(u) = (1 + ...) + B(u 2 + ...) , (4.10) 

where B is the coefficient for the second solution and the coefficient for the first solution has 
been set to 1. At the horizon the asymptotic solution satisfying the incoming wave boundary 
condition is 

= (i -u)- ir °/ 2 a + ...) . (4.H) 

As discussed in [59, 60] we find the coefficient B by integrating the two boundary solutions 
from (4.10) forward towards the horizon and by matching the linear combination of the nu- 
merical solutions $(w) = $^ um + B&2 Um to the solution (4.11) at the horizon. The imaginary 
part of the retarded correlator then is given by 



-2T 



-ImG 



R 



N 2 (nT) 4 2lmB 

IT tt> 



UJ 47T 2 

Solving (4.9) and matching the asymptotic solutions as described above, we obtain 



lim 



d 3 q q 2 

(2tt) 3 y 



-2T 



CO 



lmG R 



p2 p2 



Kq) 



N 2 T 9 67.258 . 



(4.12) 



(4.13) 



The corresponding result for the energy-momentum tensor correlator is obtained in an anal- 
ogous way but the analysis is significantly more complicated. Fortunately it has been ex- 
tensively and carefully analyzed [57]. The final result is 



lim / 



d 3 q q 2 

(2tt) 3 Y 



-2T 



ImG^ 



TT 



N 2 T 9 355.169 . 



(4.14) 
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V. SUMMARY AND DISCUSSION 

Over the duration of the lifetime of the heavy meson state the meson will loose momentum 
on average and simultaneously receive random kicks as codified by the Langevin equations 
of motion 

^ = -vop* + e(t) , (r = ^s(t - o . (5.i) 

The drag and momentum broadening rates are related by the Einstein relation 

Vd = 2TM~ ' (5 ' 2) 

with M the meson mass. Collecting the results for polarizabilities (3.34), (3.46) and force 
correlators (4.13), (4.14), and using Eq. (2.24) we obtain our principal result 



T 3 (2txT 
K ~N~ 2 \~M^ 



2 / i n \ 2 



, (67.258) + ( \- ] (355.169) 

071 J \07l J 



The finite temperature mass shift in M = 4 SYM is 



(5.3) 



= — T 



N 2 

2nT\ 3 9n ^ 



V M J 10 



Comparing these formulas with the analogous formulas in weak coupling large N QCD given 
in Eqs. (2.18) and (2.5) 

T 3 (ttT\ & 50176 
KpQCD ~ W V K~b) l215" 7r ' ( } 

and 

w -°=- T (£) 3 S- < 5 ' 6 > 

we see that the meson mass M plays the role of the inverse Bohr radius A# = (m q /2)a s CF 
in the strong coupling dipole effective Lagrangian. This is as expected for relativistic bound 
states. Below we will compare the values of the ratio k/ (5M 2 ) at strong and weak coupling. 

Perhaps the theoretically most important aspect of this work is that we have deduced 
a drag coefficient which is suppressed by N 2 in the large N limit. On the field theory 
side, this was achieved by calculating the mass shift of a meson in an external background 
(which is finite at large N), and using this information to deduce the meson couplings to the 
stress tensor and the operator Opt. We have restricted the calculation to the heavy dipole 
limit where these are the only relevant operators. Subsequently, the fluctuations of these 
operators give rise to a net force on the meson. 

On the gravity side, the meson mass shift arises as a change in the normal vibrational 
modes of the D7 brane in the presence of an external gravitational (or dilatonic) field. 
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FIG. 3: Feynman graph leading to scattering of a heavy meson in perturbation theory. This graph 
represents the gauge contribution to correlators given in Eq. (5.7). At strong coupling the effect of 
additional scatterings is to reduce the integrated value of this correlator by almost a factor of five 
in the appropriate kinematic regime. 



Although it is not manifest in the usual black hole finite temperature AdS/CFT setup, 
the gravitational field and dilatonic fields are continually fluctuating. This is encoded by 
the fluctuation dissipation theorem in the field theory which does emerge in a Kruskal 
formalism of the gauge gravity duality. (The fluctuation dissipation theorem was used to 
relate Eq. (2.10) and Eq. (2.12).) Since these fluctuating gravitational and dilatonic fields 
shift the spectrum of the D7 brane excitations, gradients in these fields give rise to a net force 
on a mesonic normal modes of the D7 brane. This discussion suggests a better understanding 
of how gravitational and dilatonic fields fluctuate in bulk, would give a straightforward 
procedure to calculate the drag of a finite mass meson. Specifically, fluctuations in the bulk 
would force motion of meson wave functions which extend into the fifth dimension. We hope 
to pursue this reasoning in the future. 

From a phenomenological perspective the current calculation was limited to very heavy 
mesons (which survive above T c ) where dipole interactions between the meson and the 
medium are dominant. It is certainly unclear if this is the relevant interaction mechanism 
above T c even for bottomonium. Furthermore, the dipole coupling between a heavy me- 
son and the medium is dominated by short distance physics which is not well modeled by 
AdS/CFT. 

However, after the gluons scatter off the heavy quark, they propagate out into plasma 
which modifies the free propagation as indicated by the Feynman graph in Fig. 3. To 
factorize this long distance dynamics from the short distance meson dynamics we form the 
ratio 

dt Vj<T 00 (y,t)T 00 (x,0)>| 

(5.7) 



((T 00 )) 2 

TFT 



m (8 M) 

which is independent of the short distance coefficients ct and cf provided the numerically 
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- N2 37.0 , (5.8) 



small dilatonic contribution is neglected 3 . It is then reasonable to use AdS/CFT to esti- 
mate to what degree strong coupling physics modifies this ratio in QCD. In the free finite 
temperature J\f = 4 theory, the result is (see Appendix A) 

(6M) 

Thus comparing the strong coupling result (5.7) with the weak-coupling result (5.8), we 
conclude that strong coupling effects actually reduce the scattering rate relative to the mass 
shift. Roughly speaking, the same strong coupling physics that is responsible for the reduc- 
tion of pressure (by a factor of 3/4) relative to the Stefan-Boltzmann prediction is at work 
here. However the effect is more pronounced since the correlator in Eq. (5.7) is dominated 
by larger values of spatial momentum q. 

Given this AdS/CFT result we expect the perturbative estimates for the rate of momen- 
tum broadening to be reduced by some factor which could be as large as a factor of five. We 
will not speculate on this factor here but simply write the perturbative momentum diffusion 
rate as 

- QCD = T{5Mf ^ , (5.9) 

with a gravitionally biased opinion that the coefficient is too large. To obtain a numerical 
estimate for the T(ls) state, we take T = 340 MeV ~ 2T C , M = 9.46 GeV, and estimate 
the mass shift as 5M ~ —10 MeV based on potential model calculations which fit lattice 
data [39, 40]. Substituting into Eq. (5.9) and Eq. (5.2) for the relaxation time tr = r)~ 1 we 
find 



GeV / T \ ( SM Y , 
k qcd = 0.025 — — ^-—- r — ^ , (5.10) 



InT ^340 MeVy' ViOMeVy' 

TRS i= 250fmMM(^V . ("I) 
R 7] D V9-46GeVy \ SM J K J 

The introduction discusses some of the significant uncertainties associated with the current 
formalism and the extraction of mass shift adopted here. Nevertheless, the AdS/CFT pre- 
diction is that strong coupling effects will actually increase this perturbative estimate of the 
relaxation time by up to a factor of five. 
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3 With the dilaton contribution the coefficient is 8.95. 
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FIG. 4: Euclidean graph corresponding to the Eq. (2.14) and Eq. (A3). 

APPENDIX A: DIFFUSION RATE IN PERTURBATION THEORY 

The purpose of this appendix is to compute k/(5M) 2 in free finite temperature field 
theory in M = 4 Super Yang Mills theory. This will permit a comparison to the strongly 
interacting results. 

1. QCD Computation 

In the interest of pedagogy we will indicate in detail how the QCD computation is per- 
formed. We will work in the limit where only the coupling to the electric field is included 
i.e., cb = 0. k is given by Eq. (2.14). We will use the Matsubara formalism though the real 
time formalism is not more difficult in this case. We will work in the Coulomb gauge where 
the propogators are 

J\t J d 3 xe"^ x (A (X)A o (0)) = p , (Al) 
J\r J d 3 xe-^ (A(^)^(O)) = k ^~ 2 6ii . (A2) 

Here we follow standard thermal field theory notation = (u n ,\s) and = (r, x) with 
k — |k|. u n — 2nnT labels the Matsubara and K • X — uo n r + k • x. The Euclidean metric is 
g E u = diag(+, +, +, +) and further explanation of Euclidean conventions is given in Ref. [61] 
. Then the Euclidean correlator corresponding to Eq. (2.14) and Fig. 4 for a single color 
index 

(K) = \tY:J^- 3 P°P° popo Vjl^l + coulomb graphs (A3) 

here P = P — K and we do not write the graphs involving coulomb lines since these do not 
contribute to the imaginary part. Performing the Matsubara sum, analytically continuing 
— iK° — > uj + ie, taking the imaginary part, and finally working in the limit that u> — > 
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yields the following result for the imaginary part of the retarded correlator 

at v'< ,£ ' ( "- k> = S A, (4%^ "p( i+ "p) 2 ^(^-^)(^) 3 ^(p-p- k ) w ■ 

(A4) 

with 

|M| 2 = (£ p £p) 2 (l + cos 2 (^ P p) . (A5) 

The details of the preceding steps can be streamlined and are found in many places; see 
Ref. [61] and the text book [62] for simple explanations. Integrating over the retarded 
correlator as required by Eq. (2.14) and multiplying by iV 2 to account for the number of 
gluons yields the following result 

c| 64^ 9 

iV 2 135 ' (Abj 

This result is the expected kinetic theory result for the rate of momentum diffusion of a 
heavy meson scattering via dipole scattering. 

2. = 4 Computation 

The free M = 4 Lagrangian is written as follows: 

C = 2tr j-^F 2 + l -\ a {-ia ■ d)\ a - l -d,X^X^ , (A7) 

where "a" is a SU(4) index and "2" is a SO (6) index. Under flavor rotation, A a transforms in 
the fundamental representation of SU(4) and X 1 transforms as the fundamental representa- 
tion of SO(6). SU(4) and SO(6) are locally isomorphic. SU(4) matrices are parameterized as 
e ii3 A (T*) A with trace norma ii zat ion tr[(T 4 ) A (T 4 ) B ] = C^ab and C 4 = 1/2. Similarly, SO(6) 
matrices are written as e*^- 4 ^ 6 ^, with trace normalization Cq = 1 [63]. The normalization 
convention adopted here has been fixed so that the AdS/CFT correspondence holds at the 
level of non-renormalized two point functions at zero temperature [64]. 
The full stress tensor is written 



with 



1 ~ y 1 Jgauge ^ y 1 Jfermion ^ K 1 ^scalar ' V^ ) 

rU = 2 tr |f^F- + gT (-\ F2 ) } > (A9) 

= 2 tr I lA (V^ + ^) A + sT Q A°(-«r ■ d)\ a ^J I , (A10) 

rU = 2 tr j^X^X, + <T } . (All) 



Here d = d — d . 
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There are three graphs for the gauge, fermion, and scalar loops which make a contribution 
to the imaginary part of the retarded correlator. The full diffusion rate is 

k = (k) a + (k)a + , (A12) 

where (k)a is due to gauge bosons, (k)\ is due to fermions and (k)x is due to scalars. In 
each case the retarded correlator can be written in the form of a phase space integral times 
a matrix element squared The matrix elements are 

\M\ 2 A = [N 2 ] £ p (l + cos(# pp )) 2 (l + cos 2 (# P p)) , (A13) 
\M\ 2 X = [4iV 2 ]4£ p 4 (l + cos(£ pp )) , (A14) 
\M\\ = [6N 2 ]E 4 p {l + cos6 pP ) 2 . (A15) 

Then integrating over the phase-space we obtain the three contributions to k 

S^-^w 1 *- (A18 > 

The final result for the momentum diffusion rate when only the stress tensor coupling is 
included is 



A K ) =N2 ^fT T °- (A19) 



'N 4 \ , t0 6232tt 5 

Similarly the mass shift for the meson in the finite temperature background is 

5M = (5M) A + (5M) X + (5M) X . (A20) 
The different components of the mass shift are 

( 5M )a = f^^' ( A21 ) 

(5M) X = |,[8iV 2 ^, (A22) 

( 6M )x = W [QN2]n ~^~- (A23) 

In each case the 7r 2 T 4 /30 is the energy density of a massless single component bose gas. The 
factor in square brackets counts the number of degrees of freedom (including spin) and a 
factor of 7/8 to account for the differences between bose and fermi distributions. Putting 
these pieces together we find the total mass shift due to coupling to the tensor 

7T 2 T 4 

SM = c T . (A24) 
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Now we finally evaluate the ratio in the free theory 

k 24928 ttT 



{SMY 675 N* (A25) 

This is the weak coupling expectation for this ratio provided the dominant coupling of the 
medium to the dipole is through the stress tensor operator. It is useful to compare this 
expectation to the strong coupling results as is done in the body of the text. 
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